Analogue of the Kubo Formula for Conductivity of Spatially Inhomogeneous Systems 

and Electric Fields. 
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The average densities of currents and charges, induced by a weak electromagnetic field in spatially inho- 
mogeneous systems, are calculated at finite temperatures. The Kubo formula for a conductivity tensor is 
generalized for spatially inhomogeneous systems and fields. The contributions containing electric fields and 
derivative from fields on coordinates are allocated. The Semiconductor quantum wells, wires and dots may 
be considered as spatially inhomogeneous systems. 



I. INTRODUCTION 

In connection with an increased interest to experimen- 
tal and theoretical study of low dimensional semiconduc- 
tor objects - quantum wells, wires and dots - a construc- 
tion of the fundamental theory of interaction of electro- 
magnetic fields with spatially inhomogeneous systems be- 
comes actual. 

R. J. Kubo [1] obtained the formula for the conductiv- 
ity tensor a a p(uj), applicable in a case of spatially homo- 
geneous systems and electric fields E(i) independent of 
spatial coordinates. This formula takes into account ex- 
actly the interaction of current carriers with the medium. 
Consequently, it is a powerful tool for the solution of con- 
crete conductivity problems in solids. 

In present work we generalize the Kubo formula on a 
case of spatially inhomogeneous systems and fields. Pre- 
viously we calculate average densities of currents and 
charges induced by electromagnetic field. 

Kubo [1] used the interaction operator of current car- 
riers with the electric field as 



U K 



(1) 



where and r, are the charge and radius - vector of i 
-th particle, respectively, E(t) is the time-dependent, but 
spatially homogeneous electric field. However, it follows 
from the Maxwell equations that the time-dependent 
electric field necessarily depends and on coordinates, so 
the use of Eq. (1) is always a certain approximation, 
if the field E depends on t. In spatially inhomogeneous 
systems the dependence of fields on coordinates can be 
essential. 

Our task consists in taking into account of heterogene- 
ity of systems and in obtaining of additional terms in the 
Kubo formula containing derivatives from electric field 
on coordinates. 



The operator of interaction of an electromagnetic field 
with a system of charged particles is expressed through 
the vector A(r, t) and scalar ip(r,t) potentials (see, for 
example, [2], page 68), but it is not expressed through 
electric E(r, t) and magnetic H(r, t) fields, except of in- 
dividual cases as, for example, a constant electric field. 
Accordingly and operators of densities of currents ji(r, t) 
and charges p\ (r, t) in linear approximation on vector and 
scalar potentials are not expressed through E(r, t) and 
H(r,i). However, average values (ji(r,i)) and (ji(r,i)) 
should be expressed through fields, as the observable val- 
ues. At finite temperatures T an average value is defined 
as [1,3,4] 



(...) = 



S P {exp{-f3H)...} 
Sp{exp{-f3H)} : 



TL is the Hamiltonian without an interaction of particles 
with a weak electromagnetic field. 

The expressions for average (0|ji(r, t)|0) 
and (0|pi(r, i)|0) for T = (when average (...) passes 
in average (0| . . . |0) on the ground state |0) are obtained 
in [5]. The present work is a direct continuation of [5]. 
The designations and many results of [5] will be used 
below. 

The problem of an expression of average (ji(r, t)) and 
(pi(r,t)) through electric and magnetic fields is essential 
also because if to express average (ji(r, t)) through vector 
and scalar potentials, this average contains the contribu- 
tion — (e/mc) (/o(r))A(r, t)), where e = ej,m = rrii are 
the charge and mass of the particle, c is the light veloc- 
ity, p(r) is the operator of the charge density in a zero 
approximation on fields (see Eq. (9) below). This contri- 
bution creates complexities at solutions of some concrete 
tasks, for example, about light reflection and absorption 
by semiconductor quantum wells. These complexities 
may be avoided if to express the average (ji(r, £)) and 
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(pi(r,t)) through electric fields and their derivatives on 
coordinates. 

The question on what kind of interaction to use - con- 
taining the vector potential or electric field - was dis- 
cussed earlier in [6] with reference to a task about light 
scattering in bulk crystals. In [6] apparently for the first 
time the reception is used of transition from the operators 
Vj of particles velocities to the operators r 4 of coordinates 
according to a quantum ratio v» = (i/h)\H,Ti], due to 
which it is possible to pass from expressions containing 
vector potential to expressions containing electric fields. 
However, as we should solve other tasks about spatially 
- inhomogeneous systems, it is necessary again to come 
back to this theme. 

Considering finite temperatures we use mathematical 
receptions offered in [1,3]. We compare our results with 
conclusions of [4] devoted to construction of the quan- 
tum theory of a spatial dispersion electric and magnetic 
susceptibilities. It is supposed below that there arc no 
charges and currents on the indefinitely removed dis- 
tances and that the fields E(r,t) and H(r,i) are equal 
on times t — > — oo what corresponds to adiabatic field 
switching. 



with electromagnetic field, including interaction with a 
strong magnetic field in the basic Hamiltonian 7i 

Titotai = H + U, (5) 

where the designations are introduced 

H 4^ P ' +l,(ri - rff) ' Pi = P i -- c A(r i ), (6) 

i 

U = Ui + u 2 , 

Ui = ~J d 3 rj(v)A(r,t) + J d 3 rp(rV(r, i), 



U 2 = 



2mc 



J d 3 rp(r)A 2 (r,t) 



(7) 



The operators of densities of currents and charges are 
also introduced 



j(r)=5>(r), 



ji( r ) = 2^( r _r *) v « +v i (5(r-r i )], v, = pt/m, 



II. THE SYSTEM HAMILTONIAN AND 
OPERATORS OF CURRENT AND CHARGE 
DENSITIES. 



P( r ) = ^Pi( r )i Pi(r) = e5(r-ri). 



Let us consider a system of TV particles with the charge 
e and the mass m in any arbitrary weak electromagnetic 
field, characterized by intensities E(r, t) and H(r, t). Let 
us introduce vector A(r, t) and scalar (fi(r,t) potentials 
so that 



E(r,i) = -\ dA ^ - grad<p{r,t), 
H(r,i) = rotA(r,t). 



(3) 



The fields are assumed classical, the gauge of potentials 
A and ip is any. For completeness of the task we shall 
admit, that the system of particles may be placed in the 
constant magnetic field H which may be strong. The 
vector potential A(r) corresponds to this field, so that 
H = rotA(r). The total Hamiltonian TLtotai is as follows 



H 



total 



2m 

V(ri,...r JV ) + e^<p(r i 



(4) 



where P.; = —ifi(d/dYi) is the generalized momentum 
operator, V{tx . . .r N ) is the potential energy, including 
interaction between particles and an external potential. 
In Eq. (4) it is necessary to take into account non- 
commutativity of Pi and A(ri), A(r i} t). Let us allocate 
in Eq. (4) the energy U of the interaction of particles 



The operators ji(r) and Pi(r) are connected by the con- 
tinuity equation 

divji(T)+pi( T )=0, pi(r) = ^[H,p(r)], (8) 

which will be used below. The operator U 2 is out of 
the frameworks of linear approximation on fields and is 
omitted below. 

In the Hciscnbcrg representation the additives to the 
operators of densities of currents and charges linear on 
potentials A(r, t) and <p(r, t) are equal 



ji a (r,t) = p(r,t)A a (r,t) 

mc 



dt'W),j a [T,t)\, 



Pi(r,t)= % - J' dt'piit')^^)}, 



(9) 



where the subscript 1 means a linear approximation on 
fields, p(r,t), j(r,t) and Ui(t) are operators in the in- 
teraction representation, for example, 

p(r,t) = e im/h p{r)e- mt/h 1 

[F, Q] = FQ — QF is the commutator of the operators F 
and Q. By substituting Eq. (7) into Eq. (9) we obtain 
for U x 
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j la {r,t) = ^^^(r^) 



mc 



^ij^'f ^MMr'.iOMr',*'), (10) 

pi(r,i) = ^y d\' J' dt'[p(r,t),jp(r',t')}A (r',t') 
f d 3 r' f dt'[p(r,t),p(r',t')Mr',t'). (11) 

J J — oo 



he 



III. AVERAGE VALUES OF INDUCED 
DENSITIES OF CURRENTS AND CHARGES. 

The case T = is considered in [5]. And operators of 
Eqs. (10) and (11) are averaged on the ground state of 
system. It is shown in [7] (page 84) that at averaging it 
is necessary to use wave functions |0) of the ground state 
without taking into account the interaction U. Trans- 
forming expressions for the averaged values (0|ji a (r, t)\0) 
and (0|/9i(r, i)|0) from [5] we obtain the results 

(0|j la (r,i)|0) = (0\j la (r,t)\0) E + (0\ji a (r,t)\0) dE/dr , 
(0|pi(r, t)\0) = (0|pi(r, t)\0) E + (0| Pl (r, t)\0) BE/dr , (12) 

where subscripts E and dE/dr designate the contribu- 
tions, containing an electric field and derivative from field 
on coordinate, respectively. Passing to consideration of 
systems at finite temperatures, we replace in Eq. (12) av- 
eraging (0| . . . |0) on the ground state by averaging (. . .}, 
defined in Eq. (2). Legality of such replacement will be 
proved below at comparison with results of the authors 
[1,3,4]. Thus, in Eq. (12) we replace averaging (0| . . . |0) 
by (...). According to [5] we have 

(h< x {r,t)) E = % -J d 3 r' y* dt'([jjr,t),d p (r\t')}) 



x Ep(r',t'), 



(13) 



(Mr, t)) E ^^J d 3 r' J dt'([p(r, t), dp(r', t')]) 

xE (r',t'), (14) 



(ji a (r,t)) dE/dr = — ( d/3 (r))^|M 
mc ar a 

~ J d 3 r' y* dt>([j a (r,t),Y 07 (r',t')}) 

x**£*>, (15) 



(pi(r,t))eB=--J d 3 r'J dt> 

x([p(r^),y, 7 (^0]) 9 % ( r r / , (16) 

where the designations are introduced 

d(r) = rp(r), Y^r) = r & j 7 (r), (17) 



i(r,i) = -c j dt'F,(r,t'). 

J — OO 



(18) 



Let us transform the obtained expressions so that the 
transition to the Kubo formula for spatially homoge- 
neous systems and electric field independent on coordi- 
nates would be seen clearly. We use the ratio [1,3] 

l -([F(t), Q(t')\) = y" dX(^f-F(t + ih\)), (19) 

true for any pair operators F and Q. Using Eq. (19) we 
obtain from Eq. (13) 



( ji a {r,t)) E = j d 3 r' f dt 1 [ dX 

J J-oo JO 

* { dd ^/ ] j a {r,t + ihX))E {r> ,t>). 

It is possible to show that 

ddp(r,t) dj 7 (r,t) 



dt 



= -rp- 



dr~, 



(20) 



(21) 



Substituting Eq. (21) in Eq. (20) and integrating on r' 
in parts we obtain 

(ji a (r,t)) E = f d 3 r' f dt' [ dX 

J J-oo JO 

rt r /3 



+ 



f d 3 r' f dt' f dX 

J J -oo Jo 



x (Y M (r', t') ]a (r, t + ihX)) dEp ^ ^ . (22) 

In agreement with Eq.(15) the expression for 
{jia(-, r -,t))g E /Q r consists of two parts. We do not trans- 
form first of them, and in the second we integrate on t' 
in parts and afterwards use Eq. (19). It results in 

(h a (r,t)) dE/dr = ^ c (d,(r))^^- 



-\ J dVjf dX(Y^(r')j a (v,ihX)) 



dap(r',t) 



9r 7 



f d 3 r' ( dt' f dX(Y fj7 (r',t')j a (r,t + ihX)) 

J J -oo Jo 

dEp{r>,t') 



(23) 
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Summing Eqs. (22) and (23) we see that last terms in 
the RHSs of both formulas are reduced. Total expression 
we break on two parts 



(jia(r,t)} = 0'ia(r,t)) (1) + 0'ia(r,t)) 



\(2) 



(24) 



so, that the first part contains an electric field, and the 
second contains a derivative from field on coordinates, i. 
e. 



(ji Q (r,i)}« = / dV / dt' / dX 



x(j / 3(r^^ , )Ja(r,^ + ^7iA))^; /3 (r^^ , ), 



(25) 



dV 



( jla (r,t))W = — (d (r)) 



d\{Y 07 (r')j a (r,ih\)) 



dap(r,t) 
dr a 
da f3 (r',t) 
drL ' 



(26) 



It is obvious, that the splitting Eq. (24) does not coincide 
with splitting Eq. (12) convenient only at T = 0. 
By similar way we obtain from Eqs. (14) and (16) 



(pi(r,t)) = <pi(r,i)>« + < Pl (r,i)) 



\(2) 



t ,/3 

dV / dt' / dA 

oo JO 



x{jp{r',t')p{r,t + in\))E {r',t'), 



(27) 



(28) 



( Pl {v,t))^=- l -jd\> J\\ 
x(^ 7 (r>(r,^A)>^£^. 



(29) 



By obtaining Eqs. (24) - (29) we have reached our main 
goal: We have allocated the basic contributions with sub- 
script (1) in average values of induced densities of cur- 
rents and charges and also have shown that the additional 
contributions with subscript (2) contain derivatives from, 
electric field on coordinates. The sense of splitting on ba- 
sic and additional contributions is that at solution of any 
tasks, in which the field E(r, t) is spatially inhomogc- 
neous, it is possible to estimate magnitudes of additional 
contributions and to determine: It is necessary to take 
them into account or it is possible to reject them. As 
well as in the case T = 0, obtained expressions contain 
the operators r, of particle coordinates unlike initial Eqs. 
(10) and (11) which do not contain these operators. 



IV. COMPARISON WITH RESULTS OF [4]. 

In [3,4] the expressions for values (ji(r, t)} are received 
at finite temperatures, but these expressions differ from 



deduced by us and above mentioned. Parallel we make 
similar calculations of the value (pi(r, t)), which was not 
considered in [1,3,4], but in the Maxwell equations it acts 
on the equal rights with average density of a current. Let 
us notice that in [1,3] it was considered homogeneous 
medium, in [4] - inhomogeneous. Solving the equation 
for the density matrix authors of [3] and [4] come to the 
formula, which may be obtained from Eq. (10), if in the 
RHS and in the LHS to realize averaging (...) determined 
in Eq. (2). It is obvious, that similar expression for 
(/9i(r, t)) may be obtained by averaging both sides of Eq. 
(11). 

Further the authors of [3] and [4] transform the ex- 
pression for the average induced current density in such 
a manner that the electric field E(r,i) appears in it and 
it contains the vector potential A(r,i) simultaneously. 
By doing similar procedure in initial expression for the 
average induced charge density we obtain the result of a 
kind of Eq. (27), in which the contribution (pi(r, t))^ is 
determined by Eq. (28) and the contribution (pi(r, t))^ 
is equal 

<PiM)> (2) = \f dV£" dX(j (r')p(r,ihX)) 

xA p (r',t). (30) 

For the contribution (ji Q (r, t)) (1 *> from the RHS of Eq. 
(24) in [3,4] the result of Eq. (25) is obtained, as well as 
the expression 



(j la (r,t))V=- — ( P (r))A a (r,t) 



+ -J dr 'J o dX(jp(T')j a (r,inX))A p (v',t). (31) 

The authors of [4] went further. With the help of Eq. 
(31) they have expressed the derivative from {ji a (r, t))^ 
through an electric field. Integrating this derivative on 
time we obtain 

0- lQ (r,i)) (2) --£(pW)a Q (r,i) 
+ \ j d 3 r' J dX(j (r')j a (r,iJiX))a (T',t) (32) 
and analogously 

<PiM)> (2) = \j dV J dX{j (r')p(r,ihX)) 

xa f3 {r',t). (33) 

Thus, the formulas for the additional contributions of two 
kinds are obtained: Eqs. (26) and (29), containing only 
derivatives of electric fields on coordinates, and Eqs. (32) 
and (33), containing the electric field itself. 

Let us show how it is possible to pass from Eq. (26) 
to Eq. (32), and from Eq. (29) to Eq. (33). In last term 
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from the RHS of Eq. (26) we integrate on r' y in parts. 
Then we use the equality 



: Mr) +r () {dj 1 {v)/dr 1 ) 



= j {v)-r p{v). (34) 
Then from Eq. (26) we have 

+ \ J d 3 r' J dX(j (r')j a (r,ihX))ap(r',t) 
~ J d 3 r'r' J dX{p(r')j a (T,ihX))a p (T',t). (35) 

In the last term we use Eq. (19) and integrate on r'. We 
obtain that this last term is equal 

-^(\j a (r),^2riffap(Ti,t)]) = -^-[{p{r))a a {r,t) 

i 

+ ( d/3(r) )^M]. (36) 

By substituting it in Eq. (35) we come to result of Eq. 
(32), as it was required to prove. We transform similarly 
Eq. (29) for {pi(r,t))( 2 \ The difference consists only 
that instead of Eq. (36) we have 

i 

and from Eq. (29) we pass to Eq. (33). 

Thus, comparing our results with results of [3] and 
[4], we have proved applicability of Eqs. (26) and (29) 
for the additional contributions in average values of the 
induced densities of current and charge containing only 
derivatives on coordinates from electric fields. 



V. THE ANALYSIS OF THE FORMULAS FOR 
THE ADDITIONAL CONTRIBUTION TO THE 
AVERAGE DENSITY OF A CURRENT. 

The value (ji a (r, i))*- 2 -* is determined in three forms 
of Eqs. (26), (31) and (32). In [4] some properties of 
this value are listed. Let us continue its research and 
we shall obtain the fourth expression for (ji a (r, i))*- 2 -* 
through derivatives of vector potential A(r, t) on coor- 
dinates. In the second term of Eq. (31) for jp(r') we use 
Eq. (34). It results in splitting of the second term into 
two parts: In first of them we integrate on in parts, in 
the second we use the formula 



jrWQ]) 



f 

Jo 



dX(QF(ihX)), 



(37) 



which follows from Eq. (19) at t' — t. By calculating the 
commutator we obtain 







d 3 r> I d\(Y 01 (r')j a {r,ih\)) — p d ^ . (38) 



For reception one more - fifth - expression for the ad- 
ditional contribution to the density of a current we apply 
Eq. (34) to the value j a (r,ih\) from Eq. (31), and also 
following updating of Eq. (37) 

-([F,Q]) = J d\(FQ(ih\)). (39) 



It results in 



e d 



+ 



l_d_ 

c dr^ 



( Jla (r,t))W = - — ^-{(d a (r))A (v,t)} 
J d 3 r' dX(j (r')Y ai (r,mX))A (r',t)\ . (40) 



From Eq. (40) it follows, that the integral on all space 
from (ji a (r,t)) (2S> is equal [4]. 



VI. THE TRANSITION TO EXPRESSIONS 
CONTAINING A MAGNETIC FIELD. 

Let us obtain the sixth expression for 0'i Q (r, t))( 2 > , in 
which we shall introduce a magnetic field. Let us break 
Eq. (38) on two parts 

0' 1Q M)> (2) - 0'i«(r, *)> ( - } + (ji«M)} (+) , (41) 



/■ ( +\\(±) e ia < w f dA p( r ^) , 9A a (r,t) 
~hj d3r ' Jo dA ( y /57( r ')ja(r,^A)) 



^ dAp(r>,t) ± dA 7 (r',t) S 



dr'^ ^ dr' fj 



(42) 



Since H = rotA, the value (ji Q (r,i))( ' is expressed 
through the magnetic field 

( Jla (r,i)>H,-A(H(r,i)xr) a (p(r)) 

+ Yc j dV J^(H(r',t) x v')p{ 3ti {v') 3a {v,ihX)), (43) 

and (ii a (r, t))^ may be expressed through the second 
derivatives from the vector potential on coordinates 



5 



(?ia(r,t)>< + >=- — 



2mc 



(p(r))r /3 r n 



g 2 ^ 
dr a dr 1 



+ TcJ d 3 r%r' s j o dX(j 7 (r')j a (r,ihX))-^y. (44) 

To deduce Eq. (44) from Eq. (42) we have acted as fol- 
lows: In Eq. (42) we have used a ratio Yp 7 (r') = r^j 7 (r'), 
and for j 1 (r') we have used Eq. (34) . Further in the term, 
containing dY 1 s{Y l )/dr' 5l we integrated on r' s in parts and 
in the term, containing c? 7 (r'), we used Eq. (39), inte- 
grated on r' and have calculated the commutator. 

It is possible to show that values (ji a (r,t))^ and 
{jia( r ,t))(~) separately have properties 

dw(jiM)> (±) =0, J d 3 r(j la (v,t))^ =0. (45) 

Let us show that the contribution (ji Q (r, t))^ may be 
expressed through the second derivatives from an electric 
field on coordinates. For this purpose we shall substitute 
the expression 

.MM') 



A(r, t) = a(r, t) — c 



r 

J —a 



dt' 



dr 



(46) 



for the vector potential (following from Eq. (3)) in Eq. 
(42). Then in values (ji a (r, t)}^ it is possible to allocate 
the contributions from scalar potential, to which we shall 
attribute index ip, i. e. 

0' lQ (r,t)> (±) - (.7io(r,t))L ±) + C?i«M)>W (47) 
At once we obtain that 

Vi a {T,t))P=0, (48) 

and 

d 3 r' j( dXiY^j^ihX))-^-^}. (49) 

In the second term of Eq. (49) we integrate twice in parts, 
at first on variable r 7 , then on variable r'g. Further we 
use the equation of a continuity Eq. (8) and Eq. (37). 
In a result we obtain that the second term from Eq. (49) 
is equal to the first with opposite sign and 

(7ia(r, *)><+> =0. (50) 

Taking into account Eqs. (48) and (50) we obtain from 
Eqs. (47) and (42) 



\(±) 



(j la (r,i))W = (j la (r,t)) E 

e ( dap(r,i) da a (r,t) 
2mc \ dr a drp 



~Y c J d 3 r> J dX{Y 01 {v')j la {v,inX)) 



( dap(r',t) <9a 7 (r',t) 



V dr. 



Or' 



(51) 



Taking into account the Maxwell equations rot~E = 
— (l/c)(dH/dt) and definition Eq. (18) of a(r,i) again 
we come to Eq. (43) for (j a (r, t))(~\ and we obtain also 



0i Q M)> (+) 

r/3 



e d 2 a p (v,t) 
(p(r))rpn- 



2mc 



dr a dr 1 



+ TcJ d3r 'j WM^MrSm^pr 1 - (52) 

Taking into account definition Eq. (18) we find that 
(Ji Q (r, t)}^ is expressed through the second derivatives 
from an electric field on coordinates. Thus, the sixth and 
last formula for the additional contribution (ji a (r, t))^ 
is determined by the sum of Eqs. (43) and (52). 



VII. THE EXCEPTION OF DIAGONAL 
ELEMENTS OF OPERATORS 



Eqs. (26) and (29) for the additional contributions 
in average induced densities of currents and charges, as 
against Eqs. (32) and (33), contain the operators of 
coordinates of particles. Really, the definitions Eq. (17) 
may be copied as 



d ( r ) = e^r i( 5(r - r;), 



(53) 



(54) 



But average magnitudes (ji(r, t)) and (pi(r, t)) should 
not depend on a point of readout of coordinates r,. It 
means that Eqs. (26) and (29) contain only non-diagonal 
elements of the operators r, , and the diagonal elements 
may be excluded. Let us show that it is so indeed. 

Let us transform Eq. (54) so that the operator 
stood only at the left 

Yf}l W = - ^ 5 Pi P( V ) + n ? 3iri ( r ) • (55) 

i 

Substituting Eqs. (53) and (55) in Eq. (26) we obtain 
0Ur,t)) (2) = / d 3 r> f dX{p(r')j a (r,ihX)} 



2mc 



xdiva(r',t) + — y^{r i0 5(r-r i ))- 

mc or a 

_I f d 3 r' [ dXy2(^^(r')ji a (r,ihX)) 
cj Jo 

dap(r',t) 



dr'^ 



(56) 
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First two terms in Eq. (56) we shall leave without 
changes, and in the last let us split the operator r 4 in 
two parts: 



r t = rf + rf d 



(57) 



where superscripts d and nd mean diagonal and non- 
diagonal contributions, respectively. The operator rf is 
determined through its matrix elements 



(n\rf\m) — (n|r i |m)(n|m), 



(58) 



where |n) are the eigen-functions of the Hamiltonian H. 
The commutativity property is obvious 

[H,vf]=0. 

Let us consider the contribution from the operator rf in 
the last term in Eq. (56) and designate this contribution 
as I«(r, t). We execute integration on rL f in parts and use 
a continuity ratio of Eq. (8). We obtain 

I a (T,t) = ~ J d 3 r' j\\J2{rt Pi(r')j a (r,inX)} 

x ap{r',t)). (59) 

As the operator rf commutes with the Hamiltonian H it 
is possible to rewrite Eq. (59) as 

I a (r,t) = -\J d 3 r' J P d\ 
x ^(MryUr^nxyf^hX^apir^t)). (60) 

i 

Let us substitute in Eq. (60) rfp — r^p — rfp. Then 
the value I a (r, t) will be split on two parts 



I a (r,t)=I' a (r,t)+lZ(r,t), 



(61) 



where 



dV / dX 



I' a (r,t) = - 

x^(p,(r'))j Q (r,i?iA)r J/ 3(inA))a / 3(r , ,i). (62) 

i 

I'i{r,t) = - c j d 3 r' J* dX 
x 2<Pi(r')ja(r, ihX)rfp(ihX))ap{r' \ t). (63) 

i 

Eq. (62) with the help of Eq. (37) turns in 



he 



(r)r i0 ,p(r')})a p (r',t). 



Integrating on r' and calculating the commutator we ob- 
tain the expression 

l' a (r,t) = -^Mr-r i ))^l, (65) 



dr a 



which is reduced with the second term in Eq. (56). There 
remains the value I'^r, t), which we shall transform using 
the continuity equation of Eq. (8) and integrating on 



variable r' 7 in parts. It results in 

/;>,*) = - / rf-v [ P dxj2i^(r')j a 

c J Jo 



(r, iUX) 



xrf dnX)) d -^. (66) 
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Using Eqs. (61), (65) and (66) we obtain the final ex- 
pression for the additional contribution to the average 
induced density of a current not containing of diagonal 
matrix elements of the operator r^: 

0WM)> (2) = ^/ d 3 r> j\x{p(r')j a (r,ihX)) 

xdiva(r',t) + - [ d 3 r [ dXy^{j^{r')j a (r,ihX) 
cj Jo 

x r$(ih\) 

-r^j^)ja(r,in\)) da ^ t) . (67) 



Similarly we obtain 



(p(r,t)) (2) = ^ J d 3 r' jyX{p(r')p(r,ihX))div a (r',t) 

+ - [ d 3 r' f dXY^( 3tl (r')p(r,ihX)rfp(ihX) 
cj Jo 

- r f j i7 (r')p(r,ihX))^^-. (68) 



VIII. THE CONDUCTIVITY TENSOR. 

Let us make the Fourier-transformation of an electric 
field 



£ a (k,w)= ( d 3 r ( dtE a (r,t)( 

J J — OO 



?(kr— tot) 



(69) 



(64) 



E a (r,t) = J d 3 kj^ dwE^ujy^-^ +c .c. 

(70) 

Average induced density of a current may be written 
down as 
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O'laM)) = ^4 J d3k J duj <? a p(Ku\r) 

x^fk.wje'^+c.c, (71) 

where cr a0 (k, w|r) is the conductivity tensor dependent 
on spatial coordinates (this designation is borrowed from 
[8])- 

Using Eqs. (25) and (67) for basic and additional con- 
tributions, respectively, into the average induced density 
of a current, we obtain 



a a/J (k,o;|r)=^(k,w|r)+^(k,w|r) ) (72) 



(2), 



where 



<7«(k,a;|r) = J d 3 r' dt dX(jp(r - r', -ihX) 

x^fr.tjje- 1 ^- 1 ) (73) 



a$(k,u,|r) = ^£ Ur>e-*< 



2lTlLd 



dX 



x(p(r - r', -ih\)j a {r)) + ^ J d 3 r' e -^' J dX 
x £>Vr(r - r', -ih\)j a (r)r# - rtf(-ih\) 

i 

xj^(r-r',-ih\)j a (r)). (74) 

In the case T = instead of Eq. (72) we have 

CT Q/3 ,o(k,w|r) = <7^g(k,u;|r) +<r^(k,w|r) 

+^ 7 /(k, W |r), (75) 

where 

a I a0 {^Mv)= i -j dV j~ dte-^'-^ 

i 

-r'Kp ~ r^Mr r')j a (r, t)\0), (76) 
o^(k, W |r) = ^52(0\r$6(T- ri )\0) 



x 



dV / dte- < ( kr, - wt > 
Jo 

^(0|j a (r,i)i i7 (r-r')r^ 

i 

-r$j i7 (r-r')j a (r,t)\0), 



(77) 



x(0|ba(r,i),p(r-r')]+|0), (78) 



[F, Q]+ — FQ + QF is the anti-commutator of two oper- 
ators. It is possible to pass from Eq. (75) to Eq. (72) if 
in Eqs. (76) - (78) to replace averaging (0| . . . |0) by (. . .) 
and to use Eq. (37). * 



IX. THE APPROXIMATION OF A SPATIALLY 
HOMOGENEOUS ELECTRIC FIELD. 

In some cases it is possible to neglect the contributions 
to average values of the induced densities of currents and 
charges containing derivatives from electric field on coor- 
dinates, i. e. to believe 



E(r,i)~E(t), 



(79) 



as it is made, for example, in [1] though, strictly speak- 
ing, spatially homogeneous field E may be only time- 
independent. In approximation Eq. (79) the Fourier- 
transformation is 

/oo 
dte luJt E a {t). (80) 
-oo 

Then it is possible to write down 

(j la (r,t)) h = — J dw(T a0 ((j\r)E (uj)e luJt + c.c 7 

(81) 

where the subscript h means a spatially homogeneous 
field. It is easily to see that 

cr Q/3 Hr) = <r al3 (k = 0,w|r). (82) 

Then with the help of Eqs. (72)-(74) we obtain 

/■OO r/3 

a af3 (u\v)= dt dX(J (-ihX)j a (r,t))e iut , (83) 
Jo Jo 

where 

J « = e H r la (84) 

i 

is the current operator. 

Eq. (83) is the generalization of the Kubo formula for 
the spatially inhomogeneous medium when the conduc- 
tivity tensor depends on r. 

Further we shall consider a case of a spatially homo- 
geneous medium in which any average values can not 



*In [5] the formula for the conductivity tensor is given 
at T = 0, which does not coincide with Eq. (75). This 
is because in [5] only diagonal elements (0|rj|0) are ex- 
cluded and the operator = r, — (0|r,|0) is introduced, 
distinguished from r™ d . 
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depend on coordinates r. Then the tensor CT Q/ g(o;|r) does 
not depend on r and we obtain from Eq. (83) 

1 roo rf3 

v a p{uj) = — / dt d\(J fj (-ih\)J a (t)y"\ (85) 
v o Jo Jo 

where Vq is the normalization volume. The obtained for- 
mula coincides with the Kubo result [1] if to replace in 
the RHS lo on — ui and to take into account that in [1] 
V = 1. 

X. THE CASE OF A CONSTANT MAGNETIC 
FIELD. 

Let us consider a case when an external weak electro- 
magnetic field is reduced to a constant in space and time- 
independent magnetic field H(r, t) = H and an electric 
field E(r, t) = 0. Let us remind that we have included 
the vector potential .4(r) , appropriate to a constant mag- 
netic field H, in the basic Hamiltonian Eq. (6). However, 
in the present section we believe H = 0, A(r) = and 
the field H = const is considered so weak that it is pos- 
sible to be limited by linear on a field contributions to 
the induced densities of currents and charges. Then, ac- 
cording to Eq. (25), the basic contribution to induced 
density of a current (ji(r, i))^ 1 ) = 0, since E(r, t) = 0. 
Let us choose the vector potential as 

A(r) = i(Hxr). (86) 

Then it is obvious from Eq. (44) that (ji(r, £))<+) = 
since it contains second derivatives from A(r) on coordi- 
nates. There is only contribution (ji(r, t))^ determined 
in Eq. (43). Thus, at H = const in linear approximation 
on the field we managed to express the density of the in- 
duced current through the magnetic field intensity. Now 
we shall exclude diagonal matrix elements of the opera- 
tors r, from expression for (ji(r, t)) at H = const. For 
this purpose let us take advantage of Eq. (67) in which 
vector a(r, t) can be replaced by vector potential A(r, t), 
since the initial expression of Eq. (26) can be replaced 
by Eq. (38). By substituting Eq. (86) in Eq. (67) we 
obtain 

0ia(r, t)> = £<(H x vf ) a 5(v r< )> 

i 

i 

+ YcIo dX ^ {iH X r r d Wi«( r ^A)>. (87) 

Let us notice that the vector rf d x r™ d ^ because pro- 
jections rf£ with different subscripts a do not commute 
among themselves, for example, 

(vfxvf) z = [vl d ,vl d }. (88) 



XI. CONCLUSION. 

Let us list the obtained basic results. It is shown that 
average values of densities of currents and charges in- 
duced by weak electromagnetic field at finite tempera- 
tures and spatially inhomogeneous systems are expressed 
through electric fields and their derivatives on coordi- 
nates. The contributions expressed through electric field 
were called "basic" and through derivatives - "addi- 
tional" . 

For the " additional" contributions to average values of 
induced densities of currents and charges six pairs of var- 
ious expressions are obtained. Two of these expressions 
for the density of a current coincide with the results of 
[4] . But the expressions from [4] contain electric fields or 
vector potentials, instead of derivatives from these values 
on coordinates, that complicates an estimation of value of 
the "additional" contributions. Generally speaking, in- 
tegrating on r' in parts it is possible to get rid of deriva- 
tives dEp(r', t)/dr' 7 passing to the formulas containing 
only fields instead of derivative from them. However, the 
opposite procedure - transition from a field to derivatives 
- is not possible always. In the " basic" contributions the 
fields are always kept. 

The sixth expression obtained in section VI for the 
"additional" contribution to average induced density of 
a current breaks up on two parts. First of them with 
an index (— ) is expressed only through magnetic field 
H(r, t), the second with an index (+) - through the sec- 
ond derivative from an electric field on coordinates. The 
similar result is obtained for the average induced density 
of charge. If the "additional" contributions are expressed 
through derivatives from electric fields, the operators Ti 
of coordinates of particles enter in the appropriate for- 
mulas necessary. It may seem that this result is absurd, 
since coordinate depends on a point of a beginning 
of readout. However, it appears that diagonal matrix 
elements (n\ri\n) do not enter into the average induced 
densities of currents and charges, that is shown in section 
VII, and the non-diagonal elements do not depend on a 
beginning point of readout . 

In section VIII the "basic" and "additional" contribu- 
tions are calculated in the conductivity tensor for spa- 
tially inhomogeneous systems and fields. In approxima- 
tion, when the electric field is homogeneous in space, but 
time-dependent (section IX), the " basic " contributions 
are kept only. For this case for the conductivity tensor, 
dependent on frequency w and coordinates r, the modi- 
fied Kubo formula is obtained which passes in the formula 
from [1] for spatially homogeneous systems. 

At last, in section X the expression for the average in- 
duced density of a current is obtained in a case when a 
weak electromagnetic field is reduced to a constant mag- 
netic field. 

It follows from Eq. (74) that the "additional" contri- 
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butions in conductivity contain a factor k-y/u). If a field 
E(r, t) is the plane wave extending with the light velocity 
c (at a monochromatic irradiation) or the wave package 
(at a pulse irradiation), then k ~ u> jc and additional con- 
tributions contain in comparison with basic a small factor 
v/c, where v is the speed of particles in system. However, 
this estimation is not correct always in case of spatially 
inhomogeneous systems, for example, in semiconductor 
quantum wells, wires or dots. It is possible to consider 
a field E(r, t) as an external or stimulating field only in 
the case of calculating the density of the induced current 
and charge in the lowest order on interaction of a field 
with system of charged particles. Such approximation 
is allowable in case of quantum wells under condition of 
[9,10] 7 r <C 7, where 7^(7) is the radiative (non-radiative) 
broadening of electronic excitations. 

Otherwise, when "/ r 3> 7, it is necessary to take into 
account interaction of a field with particles in all orders of 
the perturbation theory, and then E(r, t) is the genuine 
field within the low dimensional object. This field already 
cannot be presented as a superposition of plane waves for 
which k — lu/c. For example, genuine field strongly varies 
within a quantum well along an axis z, perpendicular to 
the well plane, if light is directed along an axis z, and the 
frequency w is in a resonance with one of discrete energy 
levels of the electronic system excitations in a quantum 
well [11,12]. Then the values kd ~ 1, where d is the 
quantum well width, instead of small factor v/c appears 
a the factor 



uod 2wcd' 

If a wave length A 3> d it can appear that the new fac- 
tor M is greater than v/c. In concrete cases one needs 
to estimate its value. In [11,12] it was supposed that 
M< 1 and the " additional" contributions to average in- 
duced densities of currents and charges (the case T = 
was considered ) were neglected. By substituting the re- 
ceived expressions for average densities of a currents and 
charges in the Maxwell equations, basically it is possi- 
ble to determine true fields inside and outside of low- 
dimension semiconductor objects. Thus, it is possible to 
calculate factors of reflection and absorption of light by 
these objects (see, for example, [11,12]). 

The authors are grateful M. Harkins for a critical read- 
ing of the manuscript. 
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